Wave equations on non-smooth space-times 

Gunther Hormann * 
^N) Michael Kunzinger^ 

Roland Steinbauer-t 

(N 

3 



Abstract 

We consider wave equations on Lorentzian manifolds in case of low 
regularity. We first extend the classical solution theory to prove global 
unique solvability of the Cauchy problem for distributional data and right 
hand side on smooth globally hyperbolic space-times. Then we turn to 
the case where the metric is non-smooth and present a local as well as 
a global existence and uniqueness result for a large class of Lorentzian 
manifolds with a weakly singular, locally bounded metric in Colombeau's 
algebra of generalized functions. 
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Q 1 Introduction 

O 

C — ■ In this note we are concerned with wave equations on Lorentzian manifolds, with 

a particular interest in the Cauchy problem in non-smooth situations. To achieve 
a self-contained presentation we first give a brief account on the classical solution 

. . theory in the smooth case, in particular discussing local well-posedness as well 

as global well-posedness on globally hyperbolic manifolds. We first extend this 
theory to the case of distributional data and right hand sides. Then we turn 

5J] to the case where the metric is non-smooth, that is we deal with normally 

hyperbolic operators with coefficients of low regularity. Actually the regularity 
class of the metric which we have in mind is below C 11 (i.e., the first derivative 
locally Lipschitz) — the largest class where standard differential geometric results 
such as existence and uniqueness of geodesies remain valid, and also below the 
Geroch-Traschen class of metrics — the largest class that allows for a consistent 
distributional treatment ([THII^B]). It is evident that no consistent distributional 
solution concept is available for these equations. Therefore we consider a large 
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class of weakly singular, locally bounded metrics in the setting of generalized 
global analysis and nonlinear distributional geometry ([JJJ1 Sec. 3.2]) based on 
Colombeau algebras of generalized functions ([Ml [IT]). In particular, we present 
a local existence and uniqueness result in the spirit of [17j and extend it to a 
global result on space-times with a generalized metric that allows for a suitable 
globally hyperbolic metric splitting. 

This line of research has drawn some motivation from general relativity: In 
[8] Chris Clarke suggested to replace the standard geometric definition of singu- 
larities by viewing them as obstructions to the well-posedness of the initial value 
problem for a scalar field. However, for many relevant examples as e.g. impul- 
sive gravitational waves, cosmic strings, and shell crossing singularities which 
have a metric of low regularity the Cauchy problem cannot consistently be for- 
mulated in distribution theory and one has to use a more sophisticated solution 
concept. In case of shell crossing singularities Clarke himself used a cleverly de- 
signed weak solution concept to argue for local solvability of the wave equation 
([!]). On the other hand Vickers and Wilson ([33]) used Colombeau generalized 
functions to show local well-posedness of the wave equation in conical space- 
times modelling a cosmic string. This result has been generalized to a class of 
locally bounded space-times in [17], also see [27] for the static case and [20] for 
an extension to non-scalar equations. 

This work is organized in the following way. In Section [2] we recall the clas- 
sical theory of normally hyperbolic operators on smooth Lorentzian manifolds, 
thereby essentially following a recent book by Bar et. al. ([I]). We extend their 
global existence and uniqueness result (Th. |2.7[) to the case of distributional 
data and right hand side in Section[3] In Scction|4]we recall the necessary back- 
ground from global analysis based on Colombeau generalized functions (|19|). 
We devote Section[5]to presenting a variant of the local existence and uniqueness 
theorem for the wave equation in weakly singular space-times of |17j . Finally, 
in Section [6] we extend this result to a global theorem for weakly singular space- 
times which allow for a suitable globally hyperbolic splitting of the metric. 

2 A review of wave equations on smooth 
Lorentzian manifolds 

In this section we present the solution theory for wave equations on smooth 
Lorentzian manifolds with smooth right hand side and data (cf. |3| below), 
mainly based on [TJ [31] (see also [13] ) ■ We will at once formulate the theory 
for normally hyperbolic operators on sections of a vector bundle E over some 
Lorentzian manifold M. 

We first fix some notations. Throughout we will assume (M,g) to be a 
smooth, connected, time-oriented Lorentz manifold and E — > M a smooth vector 
bundle over M. By T(M, E) or T{E) we denote the smooth sections of E — > M, 
T>(M,E) is the space of compactly supported sections, and T k (E) denotes the 
spaces of sections of finite differentiability. Given a finite-dimensional vector 



space W, the space of J-F-valued distributions in E, T>'(M, E, W) comprises the 
continuous linear maps V(M, E*) — > W. For example, given x € M, the delta- 
distribution 5 X is the E*-valued distribution in E given by 8 X : V(M, E*) -» E*, 
ip n- ip(x). 

A differential operator P : T(E) —> T(E) of second order is called normally 
hyperbolic if its principal symbol satisfies crp(£)U = —(£,:0 • id_E x (£ € T*M). 
If _B is equipped with a linear connection V then an important example of a 
normally hyperbolic operator is given by the connection d'Alembertian 

□ v : T(E) A T(T*M <g> E) A T(T*Af <g> T*M ® S) -^» r(E), 

where tr : T*M ® T*Af -»• K, tr(£ ® jj) = (£,77) is the metric trace. By the 
Weitzenbock formula, for any normally hyperbolic operator on T(E) there exists 
a unique connection V on £ and a unique B € T(Hom(E, E)) such that P = 
□ v + B. For any differential operator D : T(E) — > T(F), its formal adjoint 
D* : T(F*) -t T(E*) is uniquely characterized by 



i>(Dtp) dV = / (D*rp){y)dV 

for all ip e V(M, E) and all ip g V(M, F*). Here dV is the Lorentzian volume 
density on M. 

For notions from causality in Lorentz manifolds our main references are 
[281 [2J. In particular, for p, q <G Af, by p < q (resp. p « q) we mean that 
there exists a causal (resp. timclike) future directed curve from p to q. For 
ACflC M we write 

I+(A) := {17 e I 3p e A s.t. g >> p in fl} 

for the relative chronological future of A in f2, and analogous for J_(A). Also 
we write 

J+(A) := {q e Q I 3p e A s.t. g > p in fi}, 

for the relative causal future and analogous for J_ (A) . Here q > p means either 
p < q or p = q. A domain fl is called causal if its closure fl is contained in 
a geodesically convex domain fl' and if for all p, q £ fl the causal diamond 
J+(p) H J-{q) is a compact subset of fi'. The manifold M is said to satisfy the 
causality condition (CC) if there are no closed causal curves in M. It satisfies 
the strong causality condition (SCC) if for each p £ M and each neighborhood 
U of p there exists a neighborhood V CU such that no causal curve that starts 
and ends in V can leave U (i.e., there are no 'almost closed' causal curves in 
M). 

A Cauchy hypersurface is a subset S of M that is intersected by each inex- 
tendible timelike curve exactly once. M is called globally hyperbolic if it satisfies 
SCC and for all p, g S M the causal diamond J+ (p) fl J^} (q) is compact. The 
following is a very useful characterization of global hyperbolicity, due in its final 
form to [3] . 



Theorem 2.1. For any time-oriented Lorentzian manifold M, the following are 
equivalent: 

(i) M is globally hyperbolic. 

(ii) There exists a Cauchy hyper surf ace in M . 

(Hi) M is isometric to (ffi x S,—f3dt 2 + gt), where /3 — j3(t,x) is smooth and 
strictly positive, t i— > g t is a smooth family of Riemannian metrics on S , 
and each {t} x S is a (smooth) spacelike Cauchy-hyper surf ace in M . 

We remark that by a recent result ([!]) global hyperbolicity is also charac- 
terized by the condition that results from replacing SCC by CC in the above 
definition. 

Turning now to the problem of solving the initial value problem for a nor- 
mally hyperbolic differential operator on E — > M, we first consider the case 
where M = (V, ( , )) is a Lorentz vector space of dimension n (later on, the role 
of V will be played by a tangent space T X M to M). We denote by 7 : V -> R, 
j(X) := —(X,X) the quadratic form associated with ( , ). The analytic cen- 
terpiece of the entire construction that is to follow is provided by the so-called 
Riesz distributions: 

Definition 2.2. For a £ C with Re(a) > n, let 

R±{a){x) := ( C(a,nh(X)- ifX€J ± {0) 

1-a 2 -" 

where C(a,n) := , 2 a ™-n — tt- R+(cc) (R-(a)) is called advanced (retarded) 
Riesz distribution on V . 

R±(a) is continuous on V, and using the fact that DR±(a + 2) = R±{a) 
for Re(a) > n + 2, a 1— > R±(a) uniquely extends to a holomorphic family of 
distributions on all of C. For all a, supp(i?±(a)) C J±(0) and singsupp(_R±(a)) 
is contained in the boundary C±(0) of J±(0). Moreover, R±(0) = 5q. 

The next step in the construction is to transport the Riesz distributions 
onto the Lorentz manifold M. To this end, let fl be a normal neighborhood of 
x £ M and define the smooth function fi x : Q — > K by dV = fi x ■ (exp x )*(dz), 
where dz is the standard volume density on T X M, and exp^ is the exponential 
map at x. In normal coordinates around x, n x = y/\ detffyl. Now we set 
R±(a,x) := fi x ■ (exp x )*R±(a), i.e., 

V<^ e V(fl, C) : (i^(o, x), <p) = (R±(a), (^ • <p) o exp x ) . 

i?^:(a,a;) are called advanced (retarded) Riesz distributions on fi. The analytical 
properties of R±(a) carry over to the manifold setting, albeit in slightly more 
involved form, e.g. (setting T x := 7 o exp~ 
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UR9{a + 2 1 x) = ( - 1 -(nr T -2n) + l I R+(a,x) for a ^ 0. (1) 
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Moreover, R±(0,x) — S x , supp(i?±(a,x)) C J±(x) and singsupp(i?±(a,x)) C 
C%(z). 

For any normally hyperbolic operator P on T(E), our aim is to construct a 
fundamental solution in the following sense: 

Definition 2.3. Let P : T(E) — > T(E) be normally hyperbolic. A distribu- 
tion F e T>'(M,E,E*) such that PF = 5 X (i.e., (F,P*ip) = tp(x) for all 
(p € T>(M, E*)) is called a fundamental solution of P at x G M. F is called 
advanced (retarded) if supp(F) C J?r(x) (C. J_(x)). 

For example, in a Lorentz vector space V as above, R±(2) is an advanced 
(retarded) fundamental solution of □ at since Di?±(2) = R±(0) — Sq. In the 
manifold setting matters are more complicated, as already indicated in (fTl) . We 
first make the following formal ansatz for a fundamental solution on a normal 
neighborhood il of x: 

oo 

1l ± (x):=J2 V * R ±( 2 + 2k > x )> ( 2 ) 

k=0 

where V* £ T(fl,E ® E*), the so-called Hadamard- coefficients, are to be de- 
termined. By formal termwise differentiation one finds that in order for TZ±(x) 
to be a fundamental solution, the V™ have to satisfy the following transport 
equations: 

V grad r. V x - ( \ Dr * - " + 2k ) V " = 2kPV t 1 (* > 0) 

with V®(x) = ids x . The Hadamard coefficients are therefore uniquely deter- 
mined as the solution to this problem. 

Next, introducing convergence-generating factors into Q we obtain an ap- 
proximate fundamental solution lZ(x) in the following sense: 

3K ± € r(0 x(i,FH E) s.t. P(2)K(x) =8 X +K±(x,.). 

Here, E* ME denotes the exterior tensor product and P( 2 ) indicates that P acts 
on the second variable. We now use the K± as integral kernels to define for any 
continuous section u of E* over H: 



(K ± u)(x):= / K±(x,y)u(y)dV(y). 
Jn 

For fl a sufficiently small causal domain and any k, id + JC± : V (Cl,E*) — > 
T k (Cl, E*) is an isomorphism with bounded inverse given by the Neumann series 

oo 

(id + /C ± )- 1 =^(-/C ± F. 



Finally, for each <p £ V(Q,E*), we set Fg(.)[<p] := (id + /C ± )- 1 (7^ ± ( . )M) € 
T(E*). Then 

T(E*) -> E* 

is an advanced (retarded) fundamental solution for P at x. Indeed, 

P (2) F£(.)M = F^(.)[PV] = (id + /C ± )- 1 (^ ± (.)[P>]) 

= (id + /C ± )- 1 (P (2) ^ ± (.)M) = (id + /C±)- 1 (^ + /C ± ^)=^ 

In addition, for each ip € £>(fi,.E*), the map x' i-> F£(x')M is in T(CL,E*). 
It can be shown that the approximate fundamental solution lZ±(x) is in fact 
an asymptotic expansion of the true fundamental solution F±( . )(x) in a suit- 
able sense. Altogether, we obtain the following result on the solution of the 
inhomogeneous problem on small domains: 

Theorem 2.4. For any x <E M there exists a relatively compact causal neigh- 
borhood f2 such that the following holds: given v € D(Cl, E) and defining u± 
by 

(u±,<p):= f F2(x)[cp]-v(x)dV(x) (<p€V(Q,E)) 
Ja 

we have: u± £ T(fl,E), Pu± = v, and supp(u±) C J±(supp(i>)). 

Turning now to the global theory, the first step is to assure uniqueness of 
fundamental solutions. For this, we need certain restrictions on the causal 
structure of M: 

Theorem 2.5. Suppose that M satisfies the causality condition, that the rela- 
tion < is closed on M , and that the time separation function r : M X M —$■ R 
(see JMBI ) is finite and continuous. Let u £ T>'(M,E) be a solution of Pu = 
with future or past compact support (supp(v) n J± (p) compact for all p). Then 
u = 0. 

As an immediate corollary we obtain that under the above assumptions for 
each x £ M there is at most one fundamental solution for P at x with past 



(future) compact support. The causality conditions in Th. 2.5 are satisfied if M 



is globally hyperbolic. Under this assumption, by Th. 2.1 there exists a spacelike 
Cauchy hypersurface S in M. We denote by £ the future directed timelike unit 
normal vector field on S and consider the following Cauchy problem: 



(3) 



where / £ T>(M,E), and uq, U\ £ T>(S,E). One first notes that for vanishing 
/, mq, and u\, this problem only has the trivial solution. The analytical core of 



Pu = 


= / 


on M 


u = 


= U 


on 5 


V|U = 


= U\ 


on 5 



this result is the observation that for each ip £ T(E*) and each v £ F(E), one 
has 

ip ■ (Pv) - (PV) • v = div (W) (4) 

where the vector field VK is uniquely characterized by 

(W, X) = (V^V) • « - V> ' (Vxv) (A - € X(M )) . 

This allows to control the solution of the homogeneous equation Pu = by 
the Cauchy data on any Cauchy hypersurface. To prove existence of solutions 
one uses the above local theory to obtain solutions for /, uq and u\ supported 
in sufficiently small causal domains. More precisely, we call a causal domain 
$7 an RCCSV-domain (for relatively compact causal with small volume) if it is 
relatively compact and so small that vol(f2) • ||^±Hco(Q X m < 1- Then we have: 

Proposition 2.6. Let fi be an RCCSV-domain and suppose that f, uq and u% 
are compactly supported in tt (resp. VL n S). Then the corresponding Cauchy 
problem in fl is uniquely solvable. 

Combined with rather subtle causality arguments, this local result finally 
leads to the following main theorem on existence and uniqueness of solutions to 
§: 

Theorem 2.7. Let E be a vector bundle over a globally hyperbolic Lorentz 
manifold M with a spacelike Cauchy hypersurface S , and let P be normally 
hyperbolic on T(E). Then for each f £ T>(M,E) and each uq, u\ £ T>(S,E) 
there exists a unique solution u £ T(E) satisfying p|). Moreover, supp(u) C 
J M (supp(/) U supp(wo) U supp(ui)) and u depends continuously on (/, uo,u\). 

This result can immediately be utilized to show existence and uniqueness of 
fundamental solutions: 

Theorem 2.8. Under the assumptions of Th. \2. r /\ for each x £ M there exists 
a unique fundamental solution F + (x) (F-(x)) for P at x with past (future) 
compact support. These fundamental solutions satisfy 

(i) supp(F ± (x))CJ¥(x). 

(ii) \Jif£V{M,E*), x h-> F±(x)[(p] £T(E*) andP*(F ± {.)[tp}) = ip. 

The corresponding fundamental kernels are called Green operators: 

Theorem 2.9. Under the assumptions of Th. \2.1\ there exist unique Green 
operators G± : T>(M,E) — > T(M,E) satisfying 

(i) P o G± = i&D(M,£) 

(ii) G± o P\v(m,e) = idz>(M,E) 

(Hi) Vip £ V(M,E) : supp(G ± <p) C J| / (supp((^)) 

In fact, (G±ip){x) = F T (x)[(p\. Moreover, denoting by G± the Green operators 
forP*, 

[ (G* ± ip) • ip dV = f <p-(G T ij)dV (ip£V(M,E*),^£V(M,E)). 

JM JM 



3 Wave equations on smooth Lorentzian 
manifolds: the case of distributional data 

In this section we will be concerned with the global Cauchy problem in the 
case where the metric is still smooth but the data and right hand side are 
distributional. To keep the presentation simple we will restrict our attention to 
the wave operator. 



Based on Th. 2.1 we will assume that M = Kx S and the Lorentz metric 
on M is of the form A = —f3dt 2 + g t . We denote by □ the d'Alembertian w.r.t. 
A. The Cauchy problem we are considering can then be written as 

\M 

(5) 



□u = 


= / 


u(0,.) -- 


= "0 


V^(0, . ) = 


= U\ 



In the present setting, £ = -h&dt. To make this initial value problem meaningful 
in the distributional setting we suppose that / is smooth in the t-variable. More 
precisely, we assume that / £ C°°(R,£'(S)) n S'(R x S). It then follows from 
non-characteristic regularity (O Th. 8.3.1]) that any u £ V (M) with Du = / 
has ±(1,0) not in the wave front set WF(u)|(t, x ) for all (t,x). Thus by [121 
23.65.5], u £ C°°(M.,'D'(S)), so the initial value problem (pi) indeed makes sense 
for Uq, ui £ £'(S). 

To our knowledge, ^ has not been treated in full generality in the literature 
so far for /, uo, U\ as specified above. We therefore supply the necessary 
arguments. 

Lemma 3.1. (Uniqueness) There is at most one solution u £ C°°(K,2?'(5)) of 

Proof. We first note that □ is strictly hyperbolic with respect to the level sets 
of the map T:RxS^i, (t, x) i-> t (in the sense of [H Def. 23.2.3]). In fact, 
for the principal symbol a of □ we have 

o-\ {t , x) {T,0 = j^r 2 ~ K\x){S„0 ■ 

Thus fixing (t,x) £ M and £ £ T* S, £ 7^ 0, the polynomial p(z) :— o~\u x ){z, £) 

has the distinct real zeros ±W/3(i,x)/i^" 1 (£,£). 

Suppose now that u and u are solutions of (pM) and set w := u — u. By strict 
hyperbolicity it follows that there exists some neighborhood V of S such that 
w\ v = (see O 23.72.8]). Furthermore, by [22l Th. 23.2.9], WF(w) C Char(D) 
and is invariant under the Hamiltonian flow of a. Now Char(D) consists entirely 
of lightlikc directions. Thus the projection of any (maximal) bicharacteristic 
onto M is an inextcndible null-geodesic, hence intersects the Cauchy surface 
S. Since the wavefront set of w is empty in V and is transported along the 
bicharactcristics it therefore must be empty everywhere on M. Hence w £ 



C°°(M), and w = follows from the uniqueness part of Th. 2.7 □ 



Turning now to the problem of existence, it clearly suffices to treat the fol- 
lowing special cases of |5|: on the one hand, the homogeneous problem (/ = 0), 
which we denote by (CP1), and on the other hand the inhomogeneous problem 
with vanishing initial data, called (CP2). 



Turning first to (CP2), from Th. 2.9 it is straightforward to conclude that the 



Green operators G± : V(M) — >• C°°(M) continuously extend to operators from 
£'{M) -> V'(M) as transposed operators of G* T (see 01 Th. 4.3.10]). Thus 
given / e C°°(R,£'(S))n£'(R X S) we may set w := G+f to obtain Dw = f. 

Our task is thereby reduced to proving solvability of (CP1) since adding 
the solution of (CP1) with u = —w(0, . ) and u\ = — V*iw(0, . ) we obtain the 
desired solution of (CP2) (note that since supp(G + /) C J* f (supp/) by [34l 
(4.3.20)] both u and in are in £'(S) by 1, Cor. A.5.4]). 

To obtain a solution of (CP1) we first observe the following consequence 
of Q, (cf. [34, Th. 4.3.20]): Denote by G* := G* + - G*_ the propagator of 
the transposed operator □*. Then any smooth solution u of the homogeneous 
equation Ou = satisfies for all ip £ T>(M) 

ip-udV = f (VtG* (<p)) • u - G* {ip) ■ u x dA (6) 

M J S 

where u a = u\s, Wi = V,-u|s, and dA is the Riemannian surface element of the 
spacelike surface S (i.e., the Riemannian density w.r.t. go m our case). For the 
distributional Cauchy problem (CP1), we take ffify as a starting point and for 
given u , ui £ £'(S) define L(u ,ui) £ V(M) by 

(i(«o,«i), V>) ■= ("o, (V|G*( V ))| S ) - (ui,G*(<p)\ s ) ■ 

In case uq and u\ are the given Cauchy data u :— L(uo,Ui) will be the de- 
sired solution to (CP1). We start establishing this fact by first deriving an 
explicit formula for L(uq,Ui) well-suited to the (t, x)-splitting. This natu- 



rally has to involve the Green operators. First observe that by Th. 2.9 (ii) 



we have \3L(uq,ui) = 0, so the argument preceding Lemma 3.1 shows that in 



fact L(u QlUl ) £ C°°(R,V'(S)). Hence for ip £ D(R), ip x £ V(S) we may write 

{L(uo,ui),(po®ipi) = / ipo(t){L(uo,ui)(t),(pi)dt. (7) 

Jr 

Denote by F±(s,x) the fundamental solutions of D* at (s, x) according to 



Th.[2£] Setting F* >x := F*(s,x) - F*(s,s), we have (G*cp)(s,x) = (F* >x ,ip) 
for all (s,x) £ M and all tp £ T>(M). Furthermore, \3Ft S}X \ = ^( SjX ) — S/ S ,x) = 0) 
so F* >x £ C°°(R,T>'(S)) for each (s,x) £ M. Thus for each t £ R and each 
t/j £ T>(S) we obtain 

(L(«o,tti)(t),^> = (wo,V e -<F* x (t),V>| s =o) - (tti,<F %(*)^» • (8) 

We are now in the position to show that u indeed attains the Cauchy data. 
To this end ch oose sequences u , u i m D^S) that converge to uq resp. u\ in 



£'{S). By Th.^ for each j there exists a unique u (j) e C°°{M) such that 
Du (j) = 0, u u \0) = 4 J) , V;u (j) (0) = u^. 



In addition, again by O Th. 4.3.20] we have u (j) = L(uq\u[^). Moreover, by 
([8j| we obtain for all ip G V(S) and all* € K 

(««)(t),V.) = (Z^^X*)^) -> (i(«o,«i)(*),^) = {«(*),$ , (9) 

i.e., for all t, u(t) = limu^^(i) in V(S). In particular, u(0) = limu^(O) = 
limitQ = u , thereby verifying the first initial condition for u. 

To show that Vpm(0) = Ui we first observe that by Q we have for all 
ip £ V{S) and all t £ R 

(d t u^(t)^) = (L{u^\u^)'(t)^) 

-> (L(« ,tti)'(t),^) = <M*).V'). 



Now since for all t £ K we have Vgu(t) = l/yj/3(t, .) dtu(t) we obtain 

V =" (0) = 7IM ^ '" (0, = TPri " ma '" fa,(0) 

= lim V|U (j) (0) = lim u^ = u\ , 

thereby also verifying the second initial condition. 

Finally, we demonstrate that the support of the unique solution of (pi) sat- 



isfies the same inclusion relation as in the smooth case (Th. 2.7). To see this, 
we first note that for any u £ C°°(M) such that supp(Du) is compact we have 
the following generalization of (pf: 



cp-udV = / Du(G* + + G*_){ip)dV+ / (V*G*{<p))-u -G*((p)-u 1 dA (10) 

M JM JS 

(this follows by adapting the proof of [U Lemma 3.2.2]). Now suppose that u is 
the unique solution of (131) and pick a sequence u m £ V{M) such that u m {t) — > 
u(t) for each t £ R. For any y £ V(M) and each t £ R, supp{{G* + + G*_)(p(t, . )) 
is compact, so 



(nu m , (g; + G!» = / (o Um (t), (g; + G*_Mt, . )) dt 

(D«(t), (G; + G!>(i, ■)) «** = (□«, (G; + G!» 



Thus applying ( 10 ) to each u m and letting m — > oo we obtain 

( u ,^) = (/,(g; + gi)^) + ( Uo ,(v | g*(^))| s )-( Ui ,g*mi s ). 



From this and Th. 2.9 (iii), the claimed support properties of u follow. Summing 
up, we have proved: 

Theorem 3.2. Given f £ C°°(R,£ '(S)) H £'(R x S) and u , u x £ £'{S) there 
exists a unique solution u £ C°°(R,'D'(S)) of the Cauchy problem |5l). Moreover, 
supp(u) C J M (supp(/) U supp(wo) Usupp(ui)). 



10 



4 Generalized global analysis 

Colombeau algebras of generalized functions ([Ml EI]) are differential algebras 
which contain the vector space of distributions and display maximal consistency 
with classical analysis in the light of the Schwartz impossibility result (|30j). 
Here we review global analysis based on the special Colombeau algebra Q(M), 
for further details see QH [H [S3] and PH Sec. 3.2]. The basic idea of its 
construction is regularization of distributions by nets of smooth functions and 
the use of asymptotic estimates in terms of a regularization parameter. 

Let M be a smooth, second countable Hausdorff manifold. Set / = (0, 1] 
and denote by £(M) the subset of C ca (M) 1 consisting of all nets depending 
smoothly on the parameter e £ I. The algebra of generalized functions on M 
([H]) is defined as the quotient Q(M) := £ M (M)/Af(M) of moderate modulo 
negligible elements of £(M), where the respective notions are defined by the 
following asymptotic estimates. Here V denotes the space of all linear differential 
operators on M, 

£ M (M) ■- {(u £ ) £ e £{M) : VP CC M VP € V 3N : sup \Pu e (p)\ = 0( £ - N )} 
N{M) ■- {(u e ) e £ £{M) : VK CC M VP 6 V Vm : sup \Pu e (p)\ = 0(e m )}. 

p£K 

Elements of Q{M) are denoted by u = [(u e ) e ] = (u e ) e +Af(M). With com- 
ponentwise operations, Q(M) is a fine sheaf of differential algebras where the 
derivations are Lie derivatives with respect to smooth vector fields defined by 
Lxu := [{LxU e ) e ]i also denoted by X{u). 

There exist embeddings t of V(M) into Q(M) that are sheaf homomor- 
phisms and render C°°(M) a subalgebra of Q{M), Another, more coarse way of 
relating generalized functions in Q(M) to distributions is based on the notion of 
association: u £ Q(M) is called associated with v £ Q(M), u w v, if u 6 — v 6 — > 
in T>'(M). A distribution w £ V(M) is called the distributional shadow of u if 
u w l(w). 

The ring of constants in Q(M) is the space C of generalized numbers, which 
form the natural space of point values of Colombeau generalized functions. 
These, in turn, are uniquely characterized by their values on so-called (com- 
pactly supported) generalized points (equivalence classes of bounded nets (p s ) s 
of points, where (p e ) e ~ (<7e)e if d(p £ , q e ) — 0{e m ) for each m). 

An element u £ Q(M) is called globally bounded, if there exists a represen- 
tative (u E ) E and a C > such that |u e (x)| < C for all x G M and all e £ I. 

A similar construction is in fact possible for any locally convex space F: 
an analogous quotient construction in terms of the seminorms on F allows to 
assign a C-module Gf to P in a natural way ([E]). A particularly important 
special case is obtained for P = T(M,E), the space of smooth sections of a 
vector bundle E — » M (again with representatives that are supposed to depend 
smoothly on e). The resulting space Fg(M, E) := Gf then is a t/(M)-module, 
called the space of generalized sections of the vector bundle P. A convenient 
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algebraic description is as follows: 

Tg(M,E) = G(M) ® C oo (M) T(M,E) = L C ^ (M) (T(M,E*),G(M)). (11) 

Tg is a fine sheaf of finitely generated and projective ^/-modules. For the special 
case of generalized tensor fields of rank r,s we use the notation Q r s (M). We 
have: 

Gl(M) = Lg {M) {Ql{My,Gl{MY;Q{M)). 

Observe that this allows the insertion of generalized vector fields and one-forms 
into generalized tensors, a point of view which is essential when dealing with 
generalized metrics in the following sense. 

Definition 4.1. An element g ofG^iM) is called a generalized pseudo-Riemannian 
metric if it is symmetric (g(£,r]) = g{n,^) V£, T] G X(M)), its determinant detg 
is invertible in Q (equivalently, for each compact subset K of a chart in M there 
exists some m such that \ det(g e )ij(p)\ > e m for e small and all p € K), and it 
possesses a well-defined index v (the index of g e equals v for e small). 

Based on this definition, many notions from (pseudo-)Riemannian geometry 
can be extended to the generalized setting (cf. [IS]). In particular, any gen- 
eralized metric induces an isomorphism between generalized vector fields and 
one-forms, and there is a unique Levi-Civita connection corresponding to g. 
This provides a convenient framework for non-smooth pseudo-Riemannian ge- 
ometry and the analysis of highly singular space-times in general relativity (see 
e.g. [19, Ch. 5], [25). 

For the purposes of the present work we will also need the notion of a time- 
dependent generalized metric. To this end, let S be an n-dimensional smooth 
manifold and let pr 2 : R x S — > S denote the projection onto S. Also, let 
pr^X^S) be the corresponding pullback-bundlc. 

Definition 4.2. An element h 6 ^gipr^iT® S)) is called t- dependent generalized 
pseudo-Riemannian metric if h t possesses a well-defined index and if 

(i) (Symmetry) h t (£,r)) = h t (r),Q in Q(R x S) for all £, 77 € X(M). 

(ii) (Non-degeneracy) det(h t ) is strictly nonzero in the following sense: for 
any K CC R. and any L compact in some chart neighborhood on S there 
exists some m such that \det((h 6 ) t )ij(x)\ > e m for (t,x) £ K x L and e 
small. 

We conclude this section by the following globalization lemma. 

Lemma 4.3. Let u : I X M —$■ N be a smooth map and let (P) be a property 
attributable to values u(e,p) that is stable with respect to decreasing K and e 
in the following sense: if u(e,p) satisfies (P) for all p G K CC M and all e 
less than some £k > then for any compact set K' C K and any ek' < £k, u 
satisfies (P) on K' for all £ < £k' ■ Then there exists a smooth map it : I x M — > 
N such that (P) holds for all u(e,p) (e G /, p G M) and for each K CC M 
there exists some ek G / such that u(s,p) — u(e,p) for all (e,p) G (0,Sk] X K. 
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Proof. Let {K{)i be a compact exhaustion of M with K\ C K° +1 for all I and 
choose a smooth function 77 : M — > K with < 77 (x) < ert, for all x £ Ki\ K°_ Y 
(K := 0) (cf. e.g. QH Lem. 2.7.3]). Moreover, let v : M+ ->• [0, 1] be a smooth 
function satisfying v(t) <t for all t and 



£ < t < i 
3 

2 



"WH 1 *>3 2 



For (e,x) € Ix M let /i(e,x) := ?7(x)i/( -^y) and set u(e, x) := u(fi(e, x),x). 
Then u has all the required properties. D 



Remark 4.4. Lemma \4-3\ allows to globalize properties of Colombeau-type gen- 
eralized functions provided that representatives depend smoothly on e, as is as- 
sumed throughout this work. Indeed, if in the above situation it is additionally 
assumed that u is a representative of a Colombeau generalized function then by 
the very nature of the defining asymptotic estimates, u is itself moderate and in 
fact constitutes a representative of the same generalized function possessing the 
required property (P) globally on M . In particular, any (time- dependent) gen- 
eralized metric g possesses a representative {g e )e such that each g e is a smooth 
(time- dependent) metric globally on M , a fact that will repeatedly be used in 
what is to follow. 

5 Wave equations of non-smooth metrics I: The 
local theory 

In this section we present a local existence and uniqueness result (closely related 
to the one in [T7]) for the homogeneous wave equation of a class of generalized 
Lorentz metrics g, which will be extended to a global result in Section [6j We 
start by introducing this class of generalized weakly singular metrics. 

5.1 Weakly singular Lorentzian metrics 

Let g be a generalized Lorentzian metric on M. From now on we call the 
pair (M, g) a generalized space-time. To formulate asymptotic conditions on 
representatives (g 6 ) 6 of g let m be a background Riemannian metric on M and 
denote by || || m the norm induced on the fibers of the respective tensor bundle. 
To begin with we impose the following condition: 

(A) For all compact sets K, for all orders of derivative k g No and all fc-tuples 
of smooth vector fields rji, ... ,r]k and for any representative (g e ) e we have: 

sup||JzV..J^ 9e || m = 0{e- k ) (e->0) 

K 

su P ||jzv--^w 1 IL = °( £_fc ) ( £ ^°)- ( 12 ) 

K 
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A generalized metric with property (A) will be called a weakly singular met- 
ric. Note that here we use a somewhat different terminology as compared to 

D3I- 

Now with a view to formulating the local Cauchy problem of the wave oper- 
ator for such a metric we consider a local foliation of M given by the level sets 
of some non-singular function t € C°°(U), where [/CMis open and relatively 
compact. 

To exclude trivial cases we require the level sets E T = {q E U : t(q) = r} 
to be space-like with respect to all g e . In fact, we suppose a uniform variant 
of this condition which can also be viewed as a suitable generalization of the 
classical notion of time-orientability. Moreover, we will need a uniform bound 
on the covariant derivative of the normal form which, in particular, contains a 
condition on the second fundamental form of the level sets. More precisely, we 
demand: 

(B) Each p G M possesses a neighborhood U on which there exists a local time- 
function, that is a smooth function t with uniformly timelike differential 
alt =: a, i.e., 

g~ 1 (a,a) < —C < for some positive constant C (13) 

and one (hence any) representative g e and all small e. In addition, we 
have for all K CC U that 

sup||V e cr|| ro = 0(l) (e->0), (14) 

K 

where V is the covariant derivative of g. 

Let us denote the normal vector field to S r by £ (it is a generalized vector 
field defined via its representative £ £ G %(U), given by o — g e (£ e , •), i.e., £ e = 
grad e i). We observe that by pulling up the index, given (A), condition (B) 
is equivalent to sup x ||V e £ E || m = 0(1). Also for vector fields X,Y tangent to 
E T , we obtain V e a{X,Y) = Y(a(X)) - a(Vf,A) = 0- cr(norV^A), hence we 
obtain for the second fundamental form II e of the hypersurfaces E T 

\\II e \\ _ = 0(1) (e — > 0) uniformly on compact sets. 

Remark 5.1. Conditions (A) and (B) are given in terms of the s-asymptotics 
of the generalized metric. There is, however, the following close connection to 
the classical situation. Assume that we are given a space-time metric that is 
locally bounded but not necessarily C ' or of Geroch-Traschen class (i.e., the 
largest class that allows a consistent distributional treatment, see \lb\ \20$ and 
]3ty for the relation with the present setting). We may then embed this metric 
into the space of generalized metrics essentially by convolution with a standard 
mollifier (for details again see fS^j). From the explicit form of the embedding it 
is then clear that condition (A) holds. 

Condition (B), in adapted coordinates, demands somewhat better asymp- 
totics of the time- derivatives of the spatial part of the metric as well as the 
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spatial derivative of the (0,0) -component. This condition, in fact, is satisfied by 
several relevant examples as well. In particular, we have: 



• 



Conical space-times fall into our class since estimates (6) and (7) in 
for the embedded metric imply our condition (A), while (B) is immediate 
from the staticity of the metric. 

Impulsive pp-waves (in "Rosen form") 

ds 2 — —dudv + (1 + u + ) dx 2 + (1 — u + ) dy 2 
as well as expanding impulsive waves with line element 

ds 2 = 2dudv + 2v 2 



dz + -—Hdz 



2v 

satisfy conditions (A) and (B). Here u + denotes the kink function and 
H{z) is the Schwarzian derivative H = ti" /(2ti) - (3/i" 2 )/(4/i' 2 ) of some 
analytic function h(z) (which may be chosen arbitrarily) . For details see 
\18i Ch. 20]. In both cases the metric is continuous and it will obey con- 
ditions (A) and (B) when embedded with a standard mollifier, or - more 
generally - if we use any regularization that converges locally uniformly to 
the original metric. 

5.2 Local existence and uniqueness 

We start by formulating the local Cauchy problem for the wave operator on 
weakly singular space-times. Let p G M, and choose U to be an open and rela- 
tively compact neighborhood of p as in condition (B). Denote the corresponding 
foliation by E T = {q G U : t(q) — r} (t G [—7,7]) and suppose p G So =: E. 
In addition to the normal vector held £ and the normal covector field a whose 
asymptotics have already been discussed above we will need their corresponding 
normalized versions £ = [(£ e )e] = [(Ce/^e)e] an d o = [(f? e ) e ] = g(£, ■), where we 
have set V 2 = -&(£ e ,£ e ). 

We are interested in the initial value problem 



(15) 



where the initial data uq, U\ are supposed to be in <7(E). Note that this, in 
particular, includes the case of distributional initial data. We are interested in 
finding a local solution u G Q on U or an open subset thereof. 

A general strategy to solve PDEs in generalized functions is the following. 
First, solve the equation for fixed e in the smooth setting and form the net 
(u s ) £ of smooth solutions. This will be a candidate for a solution in Q, but 
particular care has to be taken to guarantee that the u e share a common domain 
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of definition and depend smoothly on e. In fact, as has recently been shown in 
[5] it suffices to verify continuous dependence on e. In the second step, one shows 
that the solution candidate (u e ) e is a moderate net, hence obtaining existence 
of a solution [(u e ) s ] in Q. Finally, to obtain uniqueness of solutions, one has to 
prove that changing representatives of the data leads to a solution that is still 
in the class [(u £ ) e ]. Note that this amounts to an additional stability of the 
equation with respect to negligible perturbations of the initial data. 

So, in the present situation we need a condition which provides us with the 
existence of a solution candidate: 

(C) For each p £ S there exists a neighborhood V C U and a representa- 
tive (g s ) e of the metric g on V such that V is, for each e, an RC CSV- 
neighborhood in (M, g e ) with SflFa spacelike Cauchy hypersurface for 
V. 



Indeed Prop. |2.6| now provides us with a solution candidate defined on V, 
that is a net (u e ) with O e u e — f e on V for some negligible net (/ e ) and moreover 
u e|snv = uo e , Vfu\s n v — Uie for some representatives (uq £ ), (ui e ) of the 
data. We note that continuous dependence of u on s follows readily from the 
construction steps detailed in Section[2J Observe that the only part that exceeds 
the classical condition for existence and uniqueness of solutions is a certain 
uniformity in e. Hcuristically this means that the light-cones of the metric g e 
neither vary to wildly with e nor collapse as e — > 0. In terms of regularizations 
of classical metrics which are locally bounded but not necessarily of C 11 or 
of Geroch-Traschen class, this condition will always be satisfied due to non- 
degeneracy of the classical metric. 

Now we may state the main result of this section. 

Theorem 5.2 (Local existence and uniqueness of generalized solutions). Let 
(M, g) be a generalized space-time with a weakly singular metric and assume 
that conditions (B) and (C) hold. Then, for each p £ S there exists an open 
neighborhood fi such that for all compactly supported uq: u i € G(^ H 0,), the 



initial value problem (15) has a unique solution u in Q(Cl). 



The core of the proof of Th. |5.2| consists of higher order energy estimates 
for the solution candidate whose existence is secured by condition (C). The 
energy estimates, which rely on conditions (A) and (B), will be carried out in a 
geometric setting using e-dependent energy momentum tensors and e-dependent 
Sobolev norms. These notions have been introduced in [331 and a suitable 
generalization of this method will be presented in the next section. 

5.3 Higher order energy estimates 

Let U D V be the neighborhoods of p given by conditions (B) and (C). The 
solution candidate (u e ) e is defined on V and we are going to estimate u £ on 
some suitable neighborhood of p. 

We start by introducing some notation. Let Vt be an open neighborhood 
of p with the property that Q C V, and such that the boundary of the region 
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Figure 1: Local foliation of space-time 



ttn{q eU : t(q) > 0} and the boundary of the region ttf\{q eU : t(q) < 0} is 
space-like w.r.t. all g e . Indeed, such a neighborhood exists by condition (A). We 
now concentrate on the forward- in-time part of Q (i.e., the part where t > 0) 
since the backward-in-time part can be dealt with analogously. To this end we 
set S T := S T n fl and denote by f2 T the open part of fi between E and S r . 
We denote the part of the boundary of il T with < t < t by Sq^ t , so that 
dQ T = So U S T U Sn, T (see Figure [l]). 

From now on we will adopt abstract index notation for (generalized) tensorial 
objects (see e.g. [H]). In particular, representatives of the metric g e and its 
inverse will be denoted by g £ ab and g^ b , respectively. In addition, to simplify the 
notation for tensors we are going to use capital letters to abbreviate tuples of 
indices, i.e., we will write Tj for T^ 1 ■■P r with |/| = r, \J\ = s. 

We now define a generalized Ricmannian metric by 

e := g + 2a ® <r, 



and use it in combination with the covariant derivative V of g to define "e- 
dcpcndent" Sobolcv norms and energies on U . Observe that by conditions 
(A) and (B) we have ||e e || m = 0(1) and ||V e e e || m = 0(1) on compact sets. 
We will also frequently need tensor products of e and use the notation e/j = 
with 171 = r = \J\. 



-piqi 



-p r q r 



Definition 5.3 (Sobolev norms and energies). Let Tj be a smooth tensor field 
and u a smooth function on U , <t < 7? and k,j G N . 



(i) We define the pointwise norm of Tj by \\Tj\ 



J KL^e 



JJrpKrpL 



T K T L 



and the higher order pointwise norm of u by |Vg u\ 2 := ||Vf; 



•V£,«| 
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(ii) On f2 T and S T we define the Sobolev norms 



V 'l"lll» .: = IE/ |V^(«)lV| ^ '»»< 



\n T ,e , . 






J2J |v e ">(«)| 2 /4 



^3=0 St 

respectively. Here \x — [(/i e ) e ] awd /i T = [(/x^) e ] denote the respective vol- 
ume forms on fl T and S T derived from g. Note that although in the second 
norm the integration is performed over the three-dimensional manifold S T 
only, derivatives are not confined to directions tangential to S T . 

(Hi) On Q we define the energy momentum tensors by (k > 0) 
T e ab '°(u) := -\gfu\ 

Tf<\u) := ( 9 rg b s d - \gfgf)el^ . . . ef - 1 «- 1 

x(v=v; 1 ...v; fc _ 1 u)(V5v^...v^_ l tt), 

(iv) Finally, the energy integrals are defined by 

E k ;ju) := J2 I T^(u)U» K"V T - (16) 

i=o Js * 

A straightforward calculation shows that the tensor fields T® b ' k (u) satisfy 
the dominant energy condition hence an application of Stokes' theorem yields 
the basic energy estimate (see e.g. [2TJ Sec. 4.3]) 



25* e («) < E^ 0<e (u) + ]T / (^ a T^(u) + Tf-»V^ b ) Me . (17) 

,-=o Ja - 

One key estimate in our approach is the equivalence of Sobolev norms and 



energies. Indeed using condition (A) and the estimate (13) one may derive f |17l 
Lemma 4.1]): 

Lemma 5.4 (Energy integrals and Sobolev norms). There exist constants C, C 
such that for each fc > and all e small 

c( v |MH T , £ ) 2 < £ T fe » < c"( v ||u||§ TrE ) 2 . (is) 

With this tool at hand we may derive the core estimate allowing to prove 
existence and uniqueness of solutions. 
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Proposition 5.5. Let (u £ ) be a solution candidate on V . Then, for each k > 1, 
there exist positive constants C k , C^, C k such that for all < r < 7 we have 



k-l 



E k T ^u e ) < 4; £ K) + ^( v n/ £ iin;, 1 s) 2 + ci'E 



1 



i=i 



-2(l+fe-j) 



ElAu e )dC 



+C' k " / El £ (u £ )dC. 



(19) 



Before sketching the proof of this statement, we draw the essential conclu- 
sions from it. Observe that the constant in front of the highest order term on 
the r.h.s. does not depend on e, hence we obtain, by an application of GronwalPs 
lemma: 



Corollary 5.6. Let (u e ) be a solution candidate on V. Then, for each k > 1, 
there exist positive constants C' k , C k , C' k such that for all < r < 7, 

E k T M) (2°) 

< 4 !e ( U£ ) + CU V ||/ £ ||^J 2 +^'2^- y J El e {u e )dc\e c '^ 



j=i 



C=o 



Consequently, if the initial energy (Eg e (u e )) E is a moderate (resp. negligible) 
net of real numbers, and (f e ) £ is negligible then 



sup (E k T . £ (u e )) e 

0<T<7 



moderate (resp. negligible). 



Now we sketch the proof of Proposition 5. 5 We have to estimate the right 



hand side of the basic energy estimate (17). Starting with the second term 



under the integral we use condition (B) to obtain 

\TfJ(u e )VU b \ < ||T e a ^'(u e )|U|V^ fe) |U < C\\T^{u e )\\ 



(21) 



Observe that by condition (A), the Riemannian metric e is O(l) hence || || m and 
j || ee are equivalent norms. We now may estimate ||T° b ' :/ (u e )||e e by the higher 
order pointwise norm of u. After integration this gives 



3=0 J ^r 



< CC\\u £ \\ k n ^J< 



(22) 



Now turning to the divergence term on the right hand side of (17) we start 
with orders k = 0, 1. Using the wave equation we find 



b.i. 



V e a Tf> u ( Ue ) = -u e V e u e and V e a T^(u s ) = f e V b E u 
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which after integration clearly can be estimated by the squares of v ||u e ||q 



and v ||/ e ||o e - Inserting this and (22) for k = 1 into (17) we obtain 



eUu^^eUu^ + cCwmi" r + Ci 



1 



Next we use Lemma |5.4| to estimate the Sobolev norm of u e in terms of its 
energy, i.e., 

( V Klla T , £ ) 2 = f{ V \KWs tl e?<K < C f T El e (u e )dC 
Jo Jo 

which gives the claim for k = 1 (with C'{ = 0). 

Finally, one has to estimate the divergence of the higher (i.e., k > 1) order 
energy momentum tensors. The general strategy is, of course, to rewrite terms 
containing the (fc + l) st order derivative of u e using the wave equation. This 
necessitates interchanging the oder of covariant derivatives, which introduces 
additional curvature terms. These can be estimated using condition (A). Ob- 
serve, however, that there also appear terms where the covariant derivatives 
falls on e. These terms of the form V^e^ J can be estimated thanks to condition 
(B). " □ 



We finally sketch the proof of Th. 5.2 We have already noted the existence 



of a solution candidate (u B ) on V. To prove that u e is moderate on J7 7 we 
start from moderateness of the data uq, u\. Inductively using the wave equation 
this translates into moderateness of the initial energies (Eg 6 (u e )) e . Now by 
Corollary 5.6 we obtain moderateness of the energies (E^ e (u E )) e for all < r < 



7. Finally, we use the Sobolev embedding theorem (together with the fact the 
volume is O(l) due to condition (A)) to estimate the sup-norm of u £ in terms 
of the Sobolev norms, which in turn can be bounded by the energies (E^. e {u e )) E 



due to Lemma 5.4 (for details see [171 Lemma 6.2]). So we see that moderateness 
of the energies implies moderateness of (u £ ) and we have proved existence of 
solutions. 

Uniqueness follows along the same lines replacing moderateness by negligi- 
bility. □ 

6 Wave equations of non- smooth metrics II: The 
global theory 



We extend the results of Th. |5.2| to establish existence and uniqueness of global 
generalized solutions. As in the classical situation we have to impose additional 
global conditions on the generalized Lorentzian metric to control causality prop- 
erties of space-time in the large. Thus, we begin by transferring the notion of 
global hyperbolicity to the setting of generalized space-times by appealing to a 



variant in terms of the metric splitting property stated in Th. 2.1 (iii). 
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Definition 6.1. Let g be a generalized Lorentz metric on the smooth (n + 
1)- dimensional manifold M. We say that (M,g) allows a globally hyperbolic 
metric splitting if there exists a C°° -diffeomorphism ip: M — > R X S, where S is 
an n-dimensional smooth manifold such that the following holds for the pushed 
forward generalized Lorentz metric A := tp*g on R X iS: 

(a) There is a representative (A £ ) ee j of X such that every X e is a Lorentz metric 
and each slice {to} x S with arbitrary to £ R is a (smooth, spacelike) Cauchy 
hypersurface for every X e (e € I). 

(b) We have the metric splitting of X in the form 

X = -/3dt 2 + h, 
where h £ ^gipr^iT^S)) is a t-dependent generalized Riemannian metric (in the 



sense of Def. J^.2) and j3 6 5(R X 5) is globally bounded and locally uniformly 
positive, i.e., for some (hence any) representative (/3 e ) of /3 and for every K CC 
R x S we can find a constant C > such that f3 e {x) > C holds for small e > 
and x £ K . 

(c) For every T > there exists a representative (h s ) ofh and a smooth complete 
Riemannian metric p on S which uniformly bounds h from below in the following 
sense: for all t £ [-T, T], x £ S , v £ T X S, and e £ I 

(h e ) t (v,v) > p(v,v). 

Remark 6.2. Observe that the basic splitting structure and the requirements 
on lower bounds for ft and h in the above definition display common features 
with the notion of regularly sliced space-times (f^ Ch. XII, Subsec. 11-4]) which 
provide sufficient conditions for global hyperbolicity in the smooth case (f]jj). 

Example 1. To obtain simple non-trivial examples of generalized space-times 
satisfying the conditions of Def. \6.1\ we consider Robertson- Walker space-times. 
First, we briefly recall the classical situation: Let (S,h ) be a connected Rie- 
mannian manifold, f: R— >]0, oo[ be smooth, and put Xi tx \ = —dt 2 + f(t) 2 (ho) x 
for every (t,x) £ R x S. Then the Lorentzian metric X on R x S is globally 
hyperbolic if and only if (S,ho) is complete (cf 'lj Lemma A. 5. 14])- More- 
over, if this is the case, then every slice {to} x S is a smooth, spacelike Cauchy 
hypersurface. 

We generalize the Robertson-Walker space-time by allowing as warping func- 
tion any f £ Q(M) that is globally bounded and locally uniformly positive and 
replace ho by a generalized Riemannian metric on S which is bounded below 
by some smooth complete Riemannian metric (to guarantee condition (c)). By 



Lemma 4-3 on globalization techniques we may pick representatives (/ E ) of f 
and (h Qe ) of h such that the smooth function f e is everywhere positive and ho e 
is a Riemannian metric on S for every e £ I . In addition, we assume that each 
ho e is complete. The generalized Lorentz metric X := — dt 2 + f 2 ho on R x S 
then trivially satisfies condition (b). Putting X £ :— —dt 2 + f 2 ho £ (s £ I) we 
obtain a representative of X. By completeness of the Riemann metric ho 6 , the 
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Lorentz metric X e is globally hyperbolic and every slice {to} x S is a Cauchy 
hypersurface for every e £ I . Thus condition (a) in Def. \6.1\ is also satisfied. 



From now on we consider only generalized space-times (M, g) which possess 
a globally hyperbolic metric splitting. To simplify notation we will henceforth 
suppress the diffcomorphism providing the splitting and assume that M =lxS 
and g = X with S and A as in the statement of Def. |6.1| Thus, the generalized 
space-time is represented by a family of globally hyperbolic space-times (M,g £ ) 
such that S = {0} x S is a Cauchy hypersurface for every g s (e £ I). 

Therefore we are provided with a suitable Cauchy hypersurface for the initial 
value problem for the wave operator □ corresponding to the generalized space- 
time metric g on M, i.e., the Cauchy problem 



(23) 



Here the unit normal vector field of S is given by £ = —rsdt and the initial 
data uq, u\ are assumed to belong to Q{S) and to have compact supports, e.g., 
arising by embedding distributional data from £'(S). 

Now the key strategy to establish a global version of Th. |5.2| on existence 
and uniqueness of solutions to the Cauchy probl em (p3| is as follows: From the 
classical existence and uniqueness result in Th. |2.7| for every e £ I we obtain 
a solution candidate defined on all of M, which again depends continuously on 
e. In this sense the global hyperbolic metric splitting of (M,g) replaces con- 
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dition (C) used in the proof of Th. 5.2 to produce a solution candidate. Then 



we aim at showing moderateness, thus existence of a generalized solution, as 
well as uniqueness by employing energy estimates as in the local constructions 
of Section p] Therefore it is appropriate to suppose also condition (A), i.e. 



g to be weakly singular. In the present situation we may translate (12 1 into 
corresponding asymptotic conditions on /? and h t . As for condition (B), we 
see that the existence of a suitable (in this case even global) foliation is a con- 
sequence of the globally hyperbolic metric splitting. Indeed we globally have 
g~ 1 (dt, dt) — —l//3 e < — C < for some positive C, which implies (p~3J) . The sec- 



ond asymptotic boundedness condition (14 1 in (B), i.e., sup^- ||V e di|| m = 0(1), 
now simply reads 

\\d(3 s \\ m = 0(l) (e->0) (24) 

uniformly on compact sets. As in the local setting this implies for the extrinsic 
curvature of the hypersurfaces {t} x S 

\\II E \\ = 0(1) (e — >• 0) uniformly on compact sets. 

We may now state the main result of this section. 

Theorem 6.3 (Global existence and uniqueness of generalized solutions). Let 
(M, g) be a generalized space-time with a weakly singular metric admitting a 
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globally hyperbolic metric splitting and assume that condition ( 14 ) holds. Then 



the Cauchy problem (23) has a unique solution u 6 Q(M) for all compactly 



supported uq,Ui G G(S). 

Sketch of proof : Let D e denote the wave operator corresponding to g £ . Then 
Th. |2.7| provides us with a global solution u 6 to D e = with Cauchy data uq e 
and u\ e , thereby defining a solution candidate. 

To prove existence we have to establish moderateness of the net {u £ ) s . 
Choose an exhaustive sequence of compact sets Kj (j € N) in 5. Then it 
suffices to show moderateness of (u e ) e on Lj :— [—j,j] x Kj for each j £ N. 

Fix j G N and choose p as in Def. |6.1| (c) with T — j. We cover Lj by 
finitely many lens-shaped regions as in Figure [T] that are, in turn, contained in 
coordinate neighborhoods. By the explicit construction given in [27l Sec. 3.3.4] 
(based on condition (A)) the heights of these lenses are uniformly bounded 
below on Lj by some S > due to the properties of p and /3. 

Employing the energy estimates from Section 4.3 we thereby derive moder- 
ateness estimates of u e on the strip [0, 5/2) x Kj from the moderateness of uo e 
and u\ £ . We may iterate this procedure to cover [0,j] x Kj in finitely many 
steps, and analogously for [— j,0] x Kj. 

By the same token negligibility of (u 0e ) and (ui E ) implies the corresponding 
property for (u £ ), which proves uniqueness. □ 

Finally, we remark that our methods also allow to treat the inhomogeneous 
equation as well as the inclusion of lower order terms, which have to satisfy 
certain asymptotic bounds, see |20j . 

Also, condition (A) is actually a little stronger than what is needed to prove 
moderateness of the solution candidate resp. negligibility in case of negligible 
data. Indeed, again by a result of [2U], it suffices to suppose (A) for k = 0. 
On the other hand we could use (A) to explicitly calculate e-power bounds of 
(derivatives) of the solution, which encode additional regularity information of 
our generalized solutions. 
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